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Abstract 

' Screening of Coulomb field of test charge in plasma with Bose condensate of electrically 

^5 I charged scalar field is considered. It is found that the screened potential contains several 

' different terms: one decreases as a power of distance (in contrast to the usual exponential 

p<J I Debye screening), some other oscillate with an exponentially decreasing envelope. Similar 

phenomenon exists for fermions (Friedel oscillations), but fcrmionic and bosonic systems have 
(— I I quite different features. Several limiting cases and values of the parameters are considered 

and the resulting potentials are presented. 

^ ; 1 Introduction 

^ I It is well known that an electric charge, Q, in plasma is screened according to the Debye law, so 

^ ' the long-ranged Coulomb field is transformed into the Yukawa type potential (see e.g. [BIS]): 

CNi : 

U(r)- ^ Q exp(-mgr) 
^^"^^'A^ 4^r ' 

lO ' where the Debye screening mass, mu, is expressed through the plasma temperature and chemical 

^ . potentials of the charged particles, see below eqs. (|12ll3p . Physical interpretation of this result 

I is evident: test charge polarizes plasma around, attracting opposite charge particles and thus the 

' electrostatic field drops down exponentially faster than in vacuum. Formally the Debye screening 



X 



appears from a pole at purely imaginary k in the photon propagator in plasma, {k"^ + IIoo)"^, 
where IIoo is the time-time component of the photon polarization operator. 

By an evident reason the screening effects were studied historically first in fermionic i.e. in 
electron-proton and in electron-positron plasma. For degenerate fermionic plasma another and 
quite striking screening behavior was found. Namely the screened potential drops down as a 
power of distance, 1/r^ in non-relativistic case and 1/r^ in relativistic case multiplied by an 
oscillating function, cos(fci?r) or sin(kFr), where kp is the Fermi momentum. This phenomenon 
is called Friedel oscillations Oil]. Usually it is prescribed to a sharp (non-analytic) cut-off of the 
Fermi distribution of degenerate electronic plasma at T = but maybe it is better to say that the 
effect is related to the logarithmic singularity of the photon polarization operator IIqo{lo = 0, fc). 
This type of screening is discussed in sec. [3] both in non-relativistic and relativistic cases for 
arbitrary, not necessarily zero, temperature. 

Plasma with charged bosons attracted attention much later, both for pure scalar electrody- 
namics, for a review see e.g. ref. [5], or for quark-gluon plasma ^ElIT]. Surprisingly until last 
year the impact of possible Bose condensate of charged fields on the photon polarization operator 
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was not considered. Only recently an investigation of plasma with Bose condensate of charged 
scalars was initiated [8]-[14j. It was found that in presence of Bose condensate the screened 
potential behaves similarly to that in fermionic case, i.e. the potential oscillates, exponentially 
decreasing with distance [9l [TO]. This effect, however, in contrast to Friedel oscillations, does 
not come from the logarithmic branch point singularity in IIoo but from the pole in the photon 
propagator at complex (not purely imaginary) value of k. It was shown that the polarization 
operator contains infrared singular term IIoo ~ |10i lllj which shifts the pole position from 
imaginary axis (as in Debye case) to a point with non-zero real and imaginary parts. 

At non-zero temperature the polarization operator has another infrared singular term ~ 1/k. 
This term is odd with respect to the parity transformation, k — > —k and, as a result, the potential 
acquires the term which decreases as a power of distance but does not oscillate. Moreover, the 
polarization operator has logarithmic singularity as in the fermionic case and this singularity also 
generates an oscillating potential similar to the Friedel one. It is interesting that the screened 
potential is a non-analytic function of the electric charge e. In particular in certain limit it may 
be inversely proportional to e, despite being calculated in the lowest order in . 

Another oscillatory, exponentially damped, behavior of the potential between static charges 
have been reported in the literature: it was argued [15] that in nuclear matter at high densities 
and low temperatures, the Debye pole acquires a non-zero real part and so the screened potential 
oscillates (see also ref. |16j). These Yukawa oscillations are short-ranged oscillations and fade 
away with distance faster, as compared to the Friedel oscillations. 

In this paper we further analyze the asymptotics of the screening effects arising in bosonic and 
fermionic plasma. In particular we have taken into account all the contributions, including the 
ones from the logarithmic singularities in the photon polarization tensor, and considered different 
limiting cases. 

The content of the paper is the following. In sec. [2] we reproduce our results for the photon 
polarization operator in plasma with charged Bose condensate. Fermionic Friedel oscillations in 
non-relativistic and relativistic cases both for T=0 and T 7^ are considered in sec. [3j There 
we reproduce already established results but using different techniques. In sec. H] we calculate 
screening in bosonic plasma, taking into account the contributions from the poles in the complex 
/c-plane, from the integral along the imaginary axis, and from the logarithmic branch cuts. The 
last part has never been done before. We consider several conditions, in particular bosons with or 
without condensate and eventually even in absence of fermions. Finally, in sec. [5] our conclusions 
are presented. 

2 Polarization operator of photon in medium 

We confine ourselves to the lowest order in the electromagnetic coupling, e^. The photon polar- 
ization operator, n^,y(a;,/c), in this approximation is well known, see e.g. books [2l[6]. For the 
calculation of the latter either imaginary or real time methods are used. However, the result can 
be obtained in a simpler way [TO] just including into the photon Green's function the effects of 
medium, namely, taking not only expectation value of the time ordered product of {A^{x)A,y{y)) 
over vacuum but add also the matter states with the weight equal to the particle distribution, 
fj{q), where j denotes the particle type and q is the particle momentum. The resulting expres- 
sions, found in many works - see e.g. and reference therein - are the following: 
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[fB{E,fi) + fB{E,fl)] X 
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1 {2q-k)^{2q-k):, 1 {2q + k)^{2q + k)^ 



2e' 



{q - kY 

d^q 
{2tt)3E 



+ 
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B 



{q + kf 



m 



B 



(2) 



[fF{E,^l) + fF{E,^l)] X 



qu{k + q)^, - q^kpg^^ + q^{k + q)^ qv{q - k)^ + q^kpg^^ + Qfiiq - k)u 



{k + qf 



+ 



rrii 



{k - qf 



nil 



(3) 



[uj, k] and q 

-2 I ^2 



[E, q] are four momenta of photon and charged particles hving in 



{q^ + p)^^"^ , with mB,F being either mass of charged bosons or fermions, 11^^ 



where k = 
plasma, E 

and Hj^^ are respectively the contribution to the polarization tensor from bosons and fermions, // 
and fl are chemical potentials for particles and antiparticles. Chemical potentials for bosons and 
fermions are generally unequal, moreover, the chemical equilibrium is not necessarily maintained 
and ^ 0. Though in what follows we present all the results for = 0, it is straightforward 
to generalize them to arbitrary n and fi. 

In kinetic equilibrium the distribution functions take the form: 



1 



B,F 



exp[(i?-/i)/r]±l' 



(4) 



where the signs " + " and " — " stay for fermions and bosons respectively. In the case that 
the boson chemical potential is equal to its maximum allowed value = niB (or Ji = niB) 
the formation of Bose condensate is possible and in equilibrium the Bose distribution function 
acquires an additional term describing accumulation of bosons in zero momentum mode: 



/B = C<5(3)(q) + 



exp [{E - mB)/T] ± 1 ' 



(5) 



where C is a constant parameter describing the amplitude of the condensate. 

The screening of test charge in the static case is determined by the zero frequency value 
of noo(0,/c). We assume that plasma is homogeneous and isotropic, so the polarization tensor 
depends only upon the magnitude of vector k but not on its direction. The corresponding 
expressions can be easily read from eqs. ()2|3p : 



Ug{0,k) 



dq ( 



27r2 Jo Eb 

e 



[fB{EB,HB) + fB{EB,fiB)] 



1 + 



E 



B 
kq 



In 



2q + k 



2q-k 



^ ^ {AEl - k^) 2q + k 



2kq 



2q-k 



(6) 
(7) 



In what follows we will omit the first argument in the polarization tensor, i.e. write Iloo{0,k) = 
Uooik). 

Evidently the first (condensate) term in /b gives rise to the quadratic infrared singularity 
Hoo ~ l/A;^ as found in refs. [KHfTT]. At non-zero temperature the pole singularity of the Bose 
distribution at g = leads to an additional infrared pole ~ 1/k in the polarization tensor of 
photons [To]. Thus at low values of the photon momentum IIoo can be expanded as [lOj : 
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c 



{2Tr)^ niB 



(8) 
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where the function h{T) is independent of k and has the hmiting values: 

h(T) - / ^'/^ ^^'^^^^ (9) 

~ I C(3/2)(mBr3)i/2/(2vr)3/2 (bwT) ' 

The low T limit of the function h{T) is however always sub-dominant with respect to the second 
term in eq. ([8]) which comes from the logarithmic term in eq. ([6]). 

In the expression of the photon polarization tensor written above the singularities of IIoo due 
to pinching of the integration contour by the poles of fBiEB,mB) and the logarithmic branch 
point in the integrand of eq. ([6]) are not taken into account. It will be done below in sec. H] 

The contribution of fermions into the polarization tensor is not infrared singular, so it is 
convenient to present the latter as 



n^o(fc) = no"o(o) + [Koik) - uUo)\ , (lo) 

where 

nS,(o) = ^ J^{f + fW + E'). (11) 

In the case of relativistic fermions with non-zero chemical potential, fi, the zero momentum 
limit of n^o is HZ]: 



n^o(O) = e'r- + ^]. (12) 




This expression is valid in the limit nip <^ ^ and mp <^ T, while in non-relativistic case for 
positive (^ — m) and small T we find: 



If < m, the polarization tensor is exponentially suppressed, IIoo ~ exp[— (m — /u)/r]. For the 
Debye mass we find the well known non-relativistic result: 

"iD = (14) 

Here, as above in the bosonic case, the singularities of 11^ due to logarithmic branch point in 
integral d?]) are not included. For that see the next section. 

The potential of a test charge, Q, modified by the plasma screening effects is given by the 
Fourier transform of the photon propagator in plasma: 

Q r d^ke-Kp{ikr) _ Q /""^ dkk"^ sm{kr) _ Q r°° dkke'^'' 

^"■^"(2^7 k^ + Uoo{k) ~2^Jo k^+Uoo{k) ~ 2n^ Jo k^ + Uoo{k) " ^ ' 

Usually the integrand in eq. (jlSp is an even function on k and the integration along the line 
of positive real k can be transformed into the contour integral in the upper complex k-plane. 
However, in the case of bosons with /x^ = the polarization operator contains and odd term 
m^r/2fc, eq. ([8]), and the usual contour transformation is not applicable. So we express integral 
(jlSp through the integral along imaginary upper fc-axis plus contribution of singularities in the 
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upper A;-plane. If IIoo is an even function of k and (A;^ +noo)^^ is regular on the imaginary A:-axis, 
the imaginary part of the integral along the imaginary axis vanishes. If the integrand has a pole 
at positive imaginary k = ikp), i.e. 

-kl + Uoo{ikD) = 0, (16) 

this poles contributes into the integral as iTT6{k—ik£)) and gives rise to the usual exponential Debye 
screening. If IIoo contains an odd part, the integral along imaginary k axis gives a contribution 
to the potential which decreases only as power of distance [10] . 

There may also be poles at complex k = kp, when both real and imaginary parts of kp are 
non-zero. Such poles have been found for plasma with charged Bose condensate [9l [10]. They 
produce oscillating behavior superimposed on the exponential decrease of the potential. It was 
argued [15] that complex poles also exist in plasma of strongly interacting particles (pions and 
nucleons) and in QCD plasma. 

There are also logarithmic singularities of IIoo (A;) a-t some non-zero 2m k and the integrals 
along the corresponding cuts also produce oscillations in the screened potential but the exponen- 
tial cut-off is much weaker, it is proportional to temperature and for zero T it becomes a power 
law one. For fermions this effect, called Friedel oscillations, is known for a long time [3l|4j, while 
for bosons a similar phenomenon has not been studied before. 



3 Friedel oscillations in fermionic plasma 

We consider here the Friedel oscillations in fermionic plasma. The non-relativistic case is discussed 
in ref. [3ll31[Il], both at zero and non-zero temperatures. The relativistic case was studied in [18]. 
In what follows all four cases are presented, considered in somewhat different way. 

Singularities of IIoo (fc) in the complex /c-plane appear when the singular points of the integrand 
in eq. ([7j) in the complex (/-plane pinch the contour of integration or coincide with the integration 
limit at q = 0. The usual calculation is done at zero temperature when the fermion distribution 
tends to theta-function and hence the integral over dq goes from zero to the Fermi momentum, 
qj. The singularity in 11^ appears when the branch points of the logarithm at A; = ±2q move 
to the integration limit at q = qp. In more general case of arbitrary temperature the integrand 
is a smooth function of q and integration goes up to infinity. The integrand has two kinds of 
singularities. First, there are poles in the distribution function fp which are situated at 

ql = [fx±i7TT{2n + l)f-ml, (17) 

where n runs from to infinity. 

The second type of singularities are branch points of the logarithm at 

Qb = ±k/2 . (18) 

The singularities of IIoo (A;) are situated at such kn for which q„ and qb coincide, g„, = qb and the 
poles, qn, and branch points, qb approach the integration contour in g-plane from the opposite 
sides. Since, according to the discussion in the previous section and eq. (fT5]) . we consider k in 
the first quadrant of the complex fe-plane, only the singularities with TZek > and Znik > 
contribute to the asymptotics of the potential, i.e. 

kn = 2qn = [(/X + iTTT{2n + 1))^ - ml] ^''^ . (19) 
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Symbolically the integral in the r.h.s. of eq. (jlSp can be written as a sum of three contributions: 

POO pkn~\~ioo 

Io= [idk] + 2TTiJ2[Res]+J2 ^> (20) 

where the first integral goes along the positive imaginary axis in /c-plane, the second one is the 
sum of the residues of the poles on the integrand (if the poles are on the imaginary axis, only a 
half of the residue is to be taken), and the third term is the integral of the discontinuity over the 
branch line of the logarithmic singularity of IloQ{k). The integration contour in complex fe- plane 
is schematically depicted in Fig. [H where only one pole and one branch-cut are included. 





o 


> 


\k_ 







Figure 1: Contour of integration in complex k-plane. 



Before calculating the singular part of IIoo let us first note that we are interested only in 
singularities in the first quadrant in fc-plane and thus only contribution from — In |2g — /c| should 
be taken. Since the absolute value of the argument can be written as the limit of e ^ of 
\2q — k\ = [{2q — k)'^ + e^]^/'^, the logarithmic contribution into U{r) is given by 



In 



k + 2q 



2q 



In 



k-2q 



(n) 

The singular part IlgQ near kn can be determined as follows. The integral along the contour 
squeezed between q^ and qb is equal to the residue of the integrand at the pole multiplied by 27ri 
plus a regular part at k = kn- The pole term near g = g„ + z is equal to 



[ln(A; -2q + ie) + ln(A: -2q- ie)] /2 ^ - \u{k - 2q - ie)/2. (21) 



1 



EnT 



exp \{En - /i)/T] + 1 zqn' 
The residue in the pole gives the singular term in IIoo equal to: 

ie^T 



(22) 



Airk 



{4E^ - k^) \n{k - 2qn 



(23) 
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where qn is the pole position given by eq. (fT9]) and En = \fq^~Vrri^ . We have neglected here 
the contributions of antiparticles assuming that the chemical potential is sufficiently large. The 
discontinuity of IIoo at the branch line k = 2g„ + iy, where y runs from zero to infinity, is equal 
to 

^^^m — -'■■'■00 ■'■^00 — 2k ' ^ ' 

where upper index " + " or " — " indicate that the value of IIoo is taken on the right or the left 
hand side of the cut. 

The contribution of this discontinuity into the asymptotic behavior of U{r), eq. p^ . is equal 

to: 



Q idy k exp(— yr + 2iqnr) (—All, 



00 j 



r(«)- 



Here k = 2g„ + iy. For fermionic plasma we can neglect y in comparison with g„ because in the 
limit of large distances y ~ 1/r. However in the bosonic case a non- vanishing contribution comes 
from sub-dominant in y terms, see below. 

Below we consider separately, firstly, relativistic and, secondly, non-relativistic cases. In 
relativistic limit En = Qn and the factor in front of logarithm, eq. ()23p . and discontinuity ()24p 
vanish at the branch point and the discontinuity becomes purely imaginary in the leading order, 
AIIgQ^ = ie^Ty. This leads to a faster decrease of the screened potential in comparison with non- 
relativistic case, instead of and to the change of phase, sin(2//r) instead of cos(2/xr). 

In relativistic case, when m <C T but /i may be large, the poles are situated at: 

En = qn=f^± mTi2n + 1) . (26) 

Since |/cp > 4|g„p > 4(/Li^ — m\), then for sufficiently large /n, ^ > uif, and low T we can neglect 
IIoo ~ e^iJ? in the denominator in comparison with and obtain: 

= llw^-^^""' = =m(2,.).-<-«)-. (27) 

For non negligible T the dominant term is that with n = and though it decreases exponentially, 
the power of the exponent may be much smaller than the standard one, eq. ([1]) with mjy = e/z/vr, 
as follows from eq. (|12p . 

At small T the result is proportional to the temperature and thus formally vanishes at T = 0. 
However, at small T the total contributions of the branch points diverges as 1/T, so summing 
up all Un we find 

TT -^TT e^gr sin(2^r) exp(-27rrr) 

167r2r3/i3 1 - exp(-47rrT) ' ^ 

For T — > and large r we can take Qn = because the effective n's are of the order of n^ff ~ 
l/(47rrT) and riT ~ 1/r <C ^. 

For very small T such that rT <C 1 we obtain: 

e^Q sin(2^r) 

^-* = ^^T^' (29) 
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in agreement with ref. [18]. However, if rT > 1, then, as we mentioned above, the screened 
potential decays exponentially similar to normal Debye screening with an important difference 
that the screening mass does not contain the electromagnetic coupling, e. On the other hand, the 
magnitude of the screened potential is proportional to e^. So formally for e = the oscillating 
potential vanishes, while the Debye one tends to the vacuum Coulomb expression. 
The ratio of the main term in the potential at T 7^ to that at T = is equal to 



U{r,T) 
U{r,T = 0) 



1 - e- 



-A-KrT 



(30) 



It is always smaller than unity, i.e. the screening is weakest at T = 0. 

Let us turn now to non relativistic limit, when mp S> T, ^ — mp <C rap, and for simplicity 
= fj, — mp ^ T. The calculations go along the same lines with evident modifications. The 
poles of the distribution function / are located at 



(/x^ - mjp) + 2i7r/ir(l + 2n) 



1/2 



y/2mpfl 



1 + 



i7r/iT(l + 2n) 



m. 



(31) 



The logarithmic singular part of IIoo corresponding to this pole is given by the same eq. (j23p 
and the discontinuity on the cut is given by eq. (j24p . An essential difference now is that the 
discontinuity does not vanish near the branch point, (4£'^ — 4(7^) = Amj^ ^ 0: 



AHoo ^ e^Tmj,/k. 

Thus the contribution of the n-th pole into the screened potential is equal to: 



Un{r) 



idy exp(2fg„r — yr) 



A:2 + nS)+(fc) 



A:2 + nS)-(A:) 



(32) 



(33) 



Here, as in the relativistic case above, k = 2g„ + iy. Neglecting k'^ in comparison with Hoo, see 
eq. (fT2|) and discussion below eq. ([26]) . we obtain: 



Unir) = „ , „ Im 

^ ' 167r2g4r-2 



2iqnT 



Qe^T 

g4^2j,2^2 



cos{2\/ 2mpjlr) exp 



-27r(2n+ 1) 



\J2mpjl 



(34) 



If temperature is not extremely small, the term with n = gives the slowest decreasing part of 
the potential, but for T ^ we need to take into account the whole sum Ucut{r) = Un{r): 



Ucutir) 



64vrV2/i2 ^os[2^2mp^^r) / ^ 

' p ' 1 — ex.p i—2-Krl y' 2nip/ fi 



exp (^—TTrTy^2mp/fi 



(35) 



Asymptotically for large r but 2iTrT 2m / jl < 1 the potential tends to 

Qe^mp cos{2qpr) 



Ucutir) 



(36) 



where qp = y/2jlmp. The result agrees with that presented in ref. [18j. The potential in eq. 
(|35p is plotted in fig. [2] as a function of distance r and temperature T for mp = 0.5 MeV 
and Hp = 0.55 MeV. Temperatures vary from 10~^ MeV and 10"^ MeV, which corresponds 
to (1.16 • 10*^ - 1.16 • 10*^) K. Distances vary from 1 MeV ^ to 100 MeV ^, corresponding to 



Figure 2: Friedel oscillations for massive fermions - see eq. 135\) - with mp = 0.5 
MeV, Hp = 0.55 MeV. Temperatures are in MeV and distances in MeV~^ . The 
exponential damping at large distance and/or temperature, as well as oscillations as 
a function of the distance r, can he seen. 



(2 • 10-1^ - 2 • 10"^) cm. The main features for the plot in the relativistic case is similar to the 
non-relativistic one. 

Note in conclusion that above we have neglected IIoo in comparison with 4g^. It is justified 
for sufficiently small e^. Otherwise one has to calculate the integral more accurately taking into 
account the mild logarithmic singularity in IIoo which goes to infinity at the branch point for 
non-relativistic fermions and goes to zero for relativistic ones. 

4 Screening in bosonic plasma 

As we have already mentioned the photon polarization tensor in presence of Bose condensate is 
infrared singular, having at small k form ([8]). The terms ~ have been found in refs. |10l lllj . 
while l/A:-term, which vanishes at T = 0, has been found in ref. [lO]. Because of l//c^ term the 
pole of the photon Green's function shifts from imaginary axis in contrast to the usual Debye 
case when the pole is purely imaginary. Due to its real part the screened potential acquires an 
oscillating factor superimposed on the exponential decrease [HI [10]. The positions of poles in 
integral (fTSj) are given by the equation /c^ + IIoo (A:) = 0, which is convenient to write as: 




(37) 



where 



(38) 
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ml = (39) 
^2 = (40) 

where h{T) is defined in eq. Q and rrij^ is the fermionic Debye mass. For relativistic fermions 
it is given by eq. (jl2p and for non-relativistic ones by eq. (jl4p . If plasma is electrically neutral 
because of the mutual compensation of bosons and fermions, the chemical potential of fermions 
is expressed through the amplitude of Bose condensate and fis = rns- However, one can imagine 
the case when there are two types of charged bosons and neutrality is achieved by the opposite 
charge densities of these bosons. In such plasma the fermionic Debye mass is zero. 

In what follows we analyze different contributions to the electrostatic potential U{r) for 
different limiting values of the parameters. In sec. 14.11 we investigate further the contribution 
from the poles in integral (|15p . In sec. 14.21 we present the contribution from the imaginary axis 
which arises when the integrand in eq. (jlSp is not an even function of k. Finally in sec. 14.31 we 
calculate the contributions from integration along the branch cuts of the logarithmic terms in 
IIoo, see eq. ([6]). The integration contour is similar to that for fermions, fig{T]but the positions of 
the poles are evidently shifted, see the following subsection. 

4.1 Contribution from poles 

At low temperatures the four roots of eq. (j37p are given by: 



^1,2,3,4 = ±^ 



e^TTiQ lb y e'^rriQ — 4e^m2 



1/2 

(41) 



As is mentioned above, we are interested only in the poles in the first quadrant in the complex 
A;-plane. If e^rriQ > Ae^rn^, all the poles are purely imaginary and the Coulomb potential is 
screened exponentially, similar to the usual Debye situation. The poles on the positive imaginary 
axis are situated at 

1/2 



The contribution of these poles into the potential is 



ki^2 = ^ ( 1 ± Vl - 4m4/e2m4 ) . (42) 



47rr k\ — k. 
In the limit of small ratio m^/emQ the potential becomes: 

Q 



^(o = j- :f . (43) 



exp ( -emor ( 1 - ) ) - -g^ exp (-m^r/mo) 



2e2mQ / / e^niQ 



(44) 



Thus for a small m2 the screening, though exponential, can be much weaker than the usual Debye 
one. 

In the opposite case, e^mg < 4:e'^m2, the poles acquire real part and now only one pole 
is situated in the first quadrant. The potential oscillates around the exponentially decreasing 
envelope [9l[T0]. The result is especially simple in the limit of large m2- 

U{r)poie = ^ exp f-^Je/2m2r] cos f^Je/2m2r] . (45) 
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More interesting situation is realized at larger temperatures, when the term m\/k in the polar- 
ization operator, eq. ()37|] is non-negligible. The contribution of the poles into the asymptotics 
of the screened potential is similar to the above considered case of low T if m2 dominates in IIoo, 
but for a small m2, e.g. if C = 0, the poles are situated at A; = e^/^(-l)^/^(m|r/2)^/^. The 
potential exponentially decreases at large distances but the power of the exponent is proportional 
to temperature and at small T the decrease of U{r) may be rather weak. 

4.2 Contribution from the integral along the imaginary axis 

Because of the odd term, m\/k, in the polarization operator the imaginary part of integral (llSp 
along the imaginary axis in the complex fc-plane is non-zero and the screened potential drops as 
a power of r: 

27r2r2 7o [-(z/r)2 + e2(mg-m|r2/z2)]2^e4^6^2/^2 

The previous expression has been obtained by substituting k = iy and then z = yr. If m2 ^ 
the dominant term at large r behaves as 

However, if the temperature is not zero and the bosonic chemical potential reaches its upper 
limit, fi = niB, but the condensate is not yet formed, the term proportional to nii dominates and 
the asymptotic decrease of the potential becomes much slower: 



U{r) = - , 3 =- , (48) 



So the formation of the condensate manifests itself by a strong decrease of screening. This effect 
may be a signal of formation of Bose condensate. 

It is interesting that the screened potential is inversely proportional to the fine structure 
constant a = e'^/Air. 



4.3 Contribution from the logarithmic branch cuts 

Let us estimate now the effects of the logarithmic singularities of IIoo on the asymptotics of 
the screened potential (analogue of the Friedel oscillations). Technically the calculations are 
similar to those made in sec. [3] but the results are noticeably different. We assume here that the 
chemical potential of bosons reaches its maximum value, /i = m^. For smaller fi there is not 
much difference between bosons and non-degenerate fermions, while for ^ = itlb new phenomena 
arise, which are absent for fermions. 

The poles in the integrand of eq. ([6]), which lead to the singularities of noo(/^) in the first 
quadrant of the complex /c-plane, are situated at 

qn = {AiTmTniB)^^'^ (1 + ii^nT /niB)^^'^ ■ (49) 

Here n runs from 1 to infinity, because there is no pole at g = since the numerator of the 
integrand is proportional to q^. 

The singularities in Hoc (A;) are situated at such k where the singularities of the integrand in 
eq. dl]) pinch the integration contour, i.e. as above, at kn = 2g„. The singular part of Hoo is 
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calculated in the same way as it has been done for fermions and is equal to the residue of the 
integrand: 



°° 2-Kk \k + 2qn + ie' ' ^ ' 



where En = \Jqn + m%. 

The discontinuity of this term across the logarithmic cut is AIlQg'*'^ = e'^TE^/k. Correspond- 
ingly the contribution of this singularity into the asymptotics of U (r) is given by: 



27r2r Jo [A;2 + n^^] [k^ + nj")] 

where /c = 2g„ + iy and E'^ = q'^ + and II^q are the values of the polarization tensor on right 
and left banks of the cut. Note that at r ^ oo the effective y is small, y ~ 1/r. 

An important difference between bosonic and fermionic cases is that the position of the 
pole for fermions, eqs. (j26|31|) . does not move to zero when T — > 0, while for bosons q^ ~ 
T. Correspondingly one can neglect Hqq in comparison with k"^, while it may be an invalid 
approximation for bosons. 

Let us first consider the case of low temperatures when IIoo is dominated by the constant 
fermionic contribution, IIq^ ~ '"^D' where m?^ is given either by eq. (jl2p or (jl3p . At large r and 
non-zero T the logarithmic contribution into the screened potential is essentially given by the 
first term with n = 1: 

Ui{r) = exp (-2v/27rmBTr) cos {2^2T:mBTr^ ■ (52) 

Here we took the relativistic limit for Ilg^. The result is easy to rewrite in non-relativistic case. 
The potential in eq. (|52p is plotted in figure O The bosonic chemical potential is taken to be 
equal to its limiting value, = f^B^ and the boson mass is assumed to be the same as the 
fermion mass in fig. [2l = mp = 0.5 MeV. Such a low mass of bosons is chosen simply for 
illustration. In realistic case charged bosons are much heavier than the charged fermions, though 
it is not excluded that there exists an unknown gauge symmetry with charged bosons lighter than 
fermions. 

The temperature in fig. Ovaries from 10"'' MeV to 0.1 MeV, corresponding to (1.16- 10^- 1.16- 
10^) K, while distances vary from 1 MeV"^ to 100 MeV-\ corresponding to (2 • 10"^^ - 2 • 10"^) 
cm. 

Figure m shows the same potential but with higher mass for bosons, tub = 100 MeV, that is 
of the order of the pion mass. The fermion mass and chemical potential are taken the same as 
above. The temperature varies in the range 10"^ — 5 • 10~^ MeV or 1.16 • 10^ — 5.8 • 10^ K and 
the distance in 10"^ < r(MeV)~^ < 10 corresponding to 2 • 10~^^ < r(cm) < 2 • 10"^''. We can 
see from these figures that if we increase the boson mass, the bosonic potential fades away faster. 

In the limit of T ^ (analogous to the discussed above Friedel case) we should take the sum 
Un because all the terms are of the same order of magnitude and n^ff ~ l/(47rmB Tr^). So 
we could expect that the sum is inversely proportional to T and the potential is non-vanishing at 
r = 0, the same as in the fermionic case. However, the summation is not so simple as previously 
because we do not deal now with geometric progression, exp(— an) but with more complicated 
function, exp(— 6-y/n). Since the effective values of n are big, we can express the sum as an integral 
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100 

Figure 3: Oscillation of the electrostatic potential in presence of bosonic plasma, see 
eq. 15^) . The boson mass is equal to the fermion one in Figure\^ ms = 0.5MeV 
and the chemical potential is = "niB- Temperatures are in MeV and distances in 
MeV-^. 



U (I 




Figure 4: Oscillation of the electrostatic potential in presence of bosonic plasma - see 
eq. 15^) . The boson chemical potential is equal to its mass, fiB = ttib = lOOMeV. 
Temperatures are in MeV and distances in MeV^^. In the picture are evident 
the oscillations due to both the temperature T and the distance r as well as the 
exponential damping in both the directions. 
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and obtain, in the leading approximation IIoo = itt-'d, that the potential is proportional to the 
temperature T and hence vanishes: 



The real part of the integral /^°° dn e^'"'^^^ written above is equal to: 
exp(-2V27rrmr) 



4r 



cos(2\/ 2TrTmr) — sin(2V 2TrTmr) ) — sin(2\/ 2-KTmr) 



irmT V / 2'nmrT 



1 



that goes to the constant value —1 in the T — > limit. So the whole expression in eq. (|53|) is 
proportional to T. 

It is important to stress that the previous result is valid in the limit Tm^r^ <^ 1, which means 
that it is applicable at small distances rs <^ l/^/rriBT. On the other hand at large distances 
r and non vanishing T one should consider the expression in eq. (|52p which is similar to the 
fermionic Priedel term in eq. (|35p but has different dependence on the coupling constant e since 
it goes like e~^, while Friedel oscillations go like e^. Hence we have non analytic dependence on 
the coupling constant e in presence of bosons. Similar dependence on was found in sec. 14. 2i 

There are also differences arising from the fact that in the limit T — > the poles of the boson 
distribution function go to zero, see eq. (j49|) . while the poles of the fermion distribution function 
tend to the non vanishing value qp, see eq. p7p . Hence Priedel oscillations for fermions start from 
their maximum amplitude at T = and then exponentially decrease with temperature, while for 
bosons the effect vanishes at T = 0, then linearly increases with T and finally exponentially 
decreases. Another consequence is that the argument of the oscillating cosine function depend 
on T for bosons but not for fermions. Hence the boson potential does not oscillate at small 
temperatures. 

At high fermionic chemical potential hf and small temperature T, the boson oscillations 
typically go to at smaller distances than the fermionic ones, which are observable at distances 
r < T. On the other hand lowering the boson mass niB the exponential damping is weaker but 
at the same time oscillations fade away. 

If the condensate is formed, Hoo would be dominated by the singular term e^m^/k"^ and 
according to eq. ()5ip the contribution of n-th branch point into the screened potential becomes: 



U^(r) = . TZe 



Q Tmj 



AJkn 



(54) 



Again, at large r and non zero T the n = 1 term is dominant. It oscillates and exponentially 
decreases according to eq. (j52]) . However, the sum T^eY^^U^ vanishes as above, eq. ([53]) . 
Probably the vanishing of U^{r) at small T in the leading order is a more general feature. At 
least the sub-leading (at small T) terms in kn and in Hoo vanish as well. If we take into account 
the imaginary part of Hoo due to the logarithmic cut, the result still remains proportional to 
a power of temperature after summation. On the other hand, as we see below, in absence of 
condensate the potential not only survives at T ^ but rises as an inverse power of T. 

Let us turn now to a more interesting though probably less realistic case when fermions are 
absent in the plasma, chemical potential of bosons is maximally allowed, /i^ = itib but the 
condensate is not formed. In the standard model a neutral system has necessarily a fermionic 
component because fermions are lighter then bosons. Anyway we can imagine systems where 
the electric charge is compensated by other heavier bosons which do not condense or models 
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with extra U{1) sector and different particle content. In this situation fermions may be absent. 
Under these conditions IIoo vanishes when T — > 0. The position of the branch points of the 
logarithm kn = 2g„ also tend to zero and the screening due to logarithmic discontinuity may 
be non- vanishing at T = 0. Indeed, let us turn again to eq. ()5ip . The integral goes along 
the contour k = kn + iy and y ~ 1/r is very small. We assume that r > Xj yJTms- Thus 
k'^ ^ k^ = IGiTmTniB- Let us now estimate IIoo ai k = k^- At small temperatures, when 
= (Eb — mB)/T g^/(2mBT), IIoo can be presented as: 



„ ,,, e^'m^p.T [ dzz 



^/SitibTz + k 



^jSrriBTz — k 



(55) 



Notice in passing that if A; < ^/8mBT, then IIoo behaves as mf/k in agreement with eqs. p7|39p . 
while at large k, k > ^/SuibT, it has the following asymptotic behavior: 

^ ^/2e^mf T3/^C(3/2) 

noo(fc) « ^57^^^ , (56) 

where C(3/2) ~ 2.6. The singular part of IIoo, eq. (j50|) . at k = kn + iy is equal to 

^1/2 22^1/2^^3/2 

n[,^^ (A:„ + iy) = g^3^,^,^/ [ln{y/8VTTnmBT) + i7r/2j . (57) 

For IIqq'* the last factor is changed to (In y/8\/ imm bT — 3i7r/2). The factor in the denominator 
of the logarithm comes from \k -\- 2tQn\ — "^l^n 

I in eq. ([50|) . 

The screened potential (fSTI) at large distances, i.e. for SnTmBr'^ > 1, is dominated by n = 1. 
One can check that |noo(^i)| > {kfl, so the latter can be neglected in the denominator of eq. 
(j5ip . Keeping in mind that we will use the result below for arbitrary n for which |noo(A;i)| > |/cf |, 
we write: 



^g2ig„r 



dxe 



In {x/S^JmBT^nTr) — in ln{x /S^^mBirnTr) + Stt'^/A 



(58) 



where x = yr. For large logarithm the leading part of the integral can be approximately evaluated 
leading to the result: 



327rQ e-2V27rTmsr 



Ui(r) = ^ — ^^77 — , sin(2v/27rTmBr) . (59) 

Note that U\{r) is inversely proportional to the electric charge and formally vanishes at T — > 0, 
but remains finite if \/TmBr is not zero. 

For smaller distances, or such small temperatures that STrTm^r^ <C 1, all n up to n-max ^ 
l/{87rTmBr'^) make comparable contributions. Thus we have to sum over n. If Umax S> 1 the 
sum can be evaluated as an integral over n. Now, for large n, k^ ^ n and may be comparable to 
noo(fcn) which, according to eq. ([57|) . drops as l/^/n. noo(fcn) would be smaller by magnitude 
than kn for 

n > no ~ 10"^ (ms/T)^/^ In^/^ (VSniBTrA . (60) 



This condition makes sense if no < n^ax or r ln^/'^(v^8rnBTV^) < 5/{Tm'^)^^^. For larger r we 
return to domination of IIoo. We should check that the condition 

rln^/^iVSruBTr^) > ^/{Tm\f'^ (61) 
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does not contradict the condition of large rimax- The latter reads 



r < I/VSttTtub. (62) 

If we neglect the logarithmic factor, both conditions would be compatible for T/tub < 4 • 10^^. 
Thus both cases of dominant noo(fcn) or k"^ can be realized depending upon relation between r, 
T, and nis- 

Let us consider smaller temperatures when |noo(^n)| > l^nl- The potential in the limit of 
small TrTmBr'^ is equal to 



U (r) = Jm 



y ne ^" 2 

(in^ ( V8mij Tr ) + 37r2/4) + ir^ In^ (VSmsTr) 



(63) 



Since the sum 



Y: ne'^'^-^ « 2 / dr? ^3,4.v™^.r, ^ _ (64) 



where rj = \/n, is real in leading order in 1/ {IQirTniBr'^), a non- vanishing contribution comes 
from the imaginary part of the numerator of the integrand and we obtain for the analogue of 
Friedel oscillations in purely bosonic case: 

IjU) ~ (65) 

2e2r2m|r6 In^ (ySmiTr) ■ 

The result has some unusual features. First, the potential decreases monotonically without any 
oscillations. Second, it is inversely proportional to the temperature, so the smaller is T, the 
larger is the potential. However, the effect exists for sufficiently small r, r ^ \ / ^/lQ^fTrnB , i.e. 
if T = O.IK and tub = IGeV the distance should be bounded from above as r <C 3 • 10"*^ cm. 
Another obstacle to realization of such screening behavior is that with fixed charge asymmetry 
Bose particle should condense and the dominant term in IIoo becomes AuibC / {211)^ . In this 
conditions we arrive to potential (|54p which vanishes at T = 0. 



5 Conclusion 

We have calculated the electrostatic potential between two test charges in plasma with electrically 
charged bosons and fermions. The new part of our consideration is an inclusion of the effects of 
the Bose condensate into the screening phenomena in plasma. To this end the chemical potential 
of bosons is taken equal to the maximally allowed value, that is to the boson mass, fi = tub- 
In this case the bosonic contribution to the time-time component of the photon polarization 
operator in plasma, noo(fc), acquires an infrared singular contribution proportional to T/k even 
before formation of the condensate and l/k'^ after formation of the condensate. Such terms 
drastically change the form of the screened potential U{r). 

All the calculations have been done in the lowest order in the electromagnetic coupling con- 
stant, e. We have imposed the condition of electric neutrality of the plasma, assuming that 
bosons and fermions compensate each other charge. We have noticed however, that the screened 
potential demonstrates a very interesting and unusual behavior as a function of temperature if 
fermions are absent. Such a situation cannot be realized in realistic equilibrium plasma because 
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the lightest charged fermions (electron/positrons) are hghter than charged bosons. However, one 
can imagine a hypothetical case of a new gauge C/(l)-symmetry with lighter charged bosons. 

We started with purely fermionic plasma and reconstructed the known Friedel oscillations 
of U{r) both in relativistic and non-relativistic limits using somewhat different technique. We 
obtained an explicit expression for U (r) at non-zero temperature which, to the best of our knowl- 
edge, is absent in the literature. 

The main part of our work is dedicated to the new phenomena created by the singularities of 
noo(^) in the complex /c-plane. We reproduced the previously obtained results that due to 
term the pole of the static photon propagator acquires non-zero real part and because of that 
the screened potential oscillates with an exponentially decreasing envelope. 

At non-zero temperature and fi = niB the polarization operator obtains an odd contribution 
with respect to the transformation k — > —k and because of that the integral along the imaginary 
axis in the complex fc-plane, which determines the asymptotic behavior of U{r), becomes non- 
vanishing. It leads to monotonic power law screening U ~ 1/r^, eq. (j47p if the condensate 
has not yet been formed. After the condensate formation the screened potential behaves as 
U{r) ~ l/(e^r^T), eq. (j48p . Such a change in the screening may be a signal of the condensate 
formation. 

We have also considered an analogue of the Friedel oscillations in the bosonic case. The 
origin of the phenomenon is the same as in the fermionic case but the resulting potential is quite 
different. The Friedel oscillations can be understood as a result of pinching the integration contour 
in the complex /c-plane by the logarithmic branch point of Hoo and the poles of the bosonic (or 
fermionic) distribution functions. However, the poles of bosonic distribution moves to zero when 
temperature tends to zero, while the fermionic ones keep a finite value. This leads to completely 
different behavior of the potential as a function of temperature. The potential vanishes when T 
goes to zero for mixed bosonic and fermionic plasma. In the case that it is dominated by the first 
pole, for large r and non-zero T, it goes as in eq. (I52p and at small T the exponential screening 
is quite mild. For purely bosonic plasma the "Friedel" part of the screening is given by eqs. 
(I58p and (I65p . If Tm^r^ is not small the potential oscillates and exponentially decreases, while 
for smaller T it does not oscillates and is proportional to l/(e^T-^). The 1/e^ behavior looks 
puzzling but one should remember that it is an asymptotic result for large distances. However, if 
we take the formal limit e — > the screening would disappear together with e. Similar reasoning 
is applicable to 1/T^ behavior: this is true only for large but simultaneously sufficiently small 
distances r < l/x/lGTrm^T, when the k'^ part of the photon Green's function is sub-dominant. 

We see that the screening is quite different in different limits and it would be very interesting 
to study this rich behavior experimentally. 
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